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Abstract 

We present a classification of SU{2>)c x SU{2)l x U{1)y gauge invariant 
AL = 2 (L being lepton number) effective operators relevant for generating 
small Majorana neutrino masses. Operators of dimension up to 11 have been 
included in our analysis. This approach enables us to systematically identify 
interesting neutrino mass models. It is shown that many of the well-known 
models fall into this classification. In addition, a number of new models are 
proposed and their neutrino phenomenology is outlined. Of particular interest 
is a large class of models in which neutrinoless double beta decays arise at a 
lower order compared to the neutrino mass, making these decays accessible 
to the current round of experiments. 
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I. INTRODUCTION 



Evidence for small but nonzero neutrino masses has been mounting over the years from 
experiments on atmospheric [|I| , solar and accelerator neutrinos. All experiments point 
toward a spectrum where the neutrino masses are less than or of order 1 eV. This is the first 
experimental challenge to the Standard Model, for its structure and particle content, taken 
as a stand-alone renormalizable field theory, do not allow for nonzero neutrino masses. New 
ingredients must be introduced in order to accommodate massive neutrinos; and the race is 
on to find the underlying new physics. 

While it is possible to accommodate neutrino masses, rriu, by modifying the low-energy 
particle content of the Standard Model, the smallness of m^, will remain a puzzle in such 
attempts. Consider, for example, introducing an SU (2) l triplet scalar field which acquires 
a vacuum expectation value (VEV) and generates Majorana masses for the left-handed 
neutrinos Either the VEV of the triplet or the relevant Yukawa couplings must be 
exceedingly small in order to generate m,^ of order 1 eV. Similarly, if right-handed neutrinos 
are introduced into the low energy spectrum so that Dirac masses can be generated, the 
smallness of the relevant Yukawa couplings (~ 10^^^) would beg for an explanation. 

A more natural explanation for the smallness of m^^ is that they are generated (via 
some underlying new physics) at a scale A higher than the electroweak scale|^ and manifest 
themselves at low energies through effective higher dimensional {d > 4) operators which are 
suppressed by appropriate powers of A. One can then understand on purely dimensional 
ground why the m,^ are small, without precise knowledge of the underlying new physics. 

As an example of this effective operator description, consider the well-known seesaw 
mechanism [^]. Here the underlying theory has a modified particle content, viz., the addition 
of heavy right-handed neutrinos. Upon integrating out these heavy fields, one arrives at an 
effective theory without the right-handed neutrinos, but with a set of dimension 5 operators 



^ Typically, A corresponds to the scale at which lepton number conservation is violated. 
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IP which can generate small Majorana masses for the left-handed neutrinos. 

While the d = 5 seesaw operators are the lowest dimensional effective neutrino mass 
operators, there may be situations in which they are not the suitable ones. For instance, as 
exemplified in many of the models discussed in Sec. Ill, there may be selection rules which 
forbid their presence. Such selection rules may be necessary for models in which the lepton 
number (L) breaking scale A is sufficiently lowQthat the d = 5 operators would generate too 
large a neutrino mass (larger than 0(1) eV). In such cases, the dominant contributions to 
the neutrino mass will come from operators with dimension higher than 5, typically through 
radiative corrections. In view of the current interests in neutrino mass models, it will be 
useful to identify all potentially relevant neutrino mass operators. We list in this paper all 
such effective operators and estimate the size of the neutrino mass they generate. In so doing 
we will reproduce many of the neutrino mass models that already exist in the literature and 
find a systematic way to arrive at new models. It should be noted that this effective operator 
approach has been widely used in other studies of possible new physics in particular, 
baryon number violation and proton decays 0,^. 

Since we are interested in operators that can lead to a mass term for the left-handed 
neutrino fields in the Standard Model, the operators must violate lepton number by two 
units, i.e., AL = 2. They contain only fields which are present in the Standard Model and 
must be SU{3)c x SU{2)l x f/(l)y invariant. They are also required to conserve baryon 
number, otherwise the nonobservation of proton decays will constrain the scale A to be 
greater than 10^^ GeV, in which case neutrino masses induced by operators with d > 5 will 
be too small to be of interest. Because limits on lepton flavor violating processes such as 
fi — i> eee and /x — >■ 67 constrain A to be larger than a few TeV, we shall focus on effective 

^For example, if L is broken due to quantum effects of gravity, A will be of order the Planck scale 
which can be as low as a few TeV in some of the scenarios that speculate the existence of large 
extra dimensions. 
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operators of dimension 11 and lower: Operators of higher dimension will likely lead to 
neutrino masses that are too small to satisfy the atmospheric neutrino data which require 



The remainder of this paper is organized as follows. Sec. II provides a classification of 
AL = 2 operators with dimension less than 12. In Sec. Ill, we present various renormaliz- 
able models that induce the AL = 2 operators of Sec. II. There we also discuss the radiative 
neutrino mass generation mechanisms in these models and outline their main phenomeno- 
logical consequences. In Sec. IV, we point out the significance of some of the operators 
to neutrino less double beta (/3/30z/) decays. We identify a number of effective operators 
which induce f3 (30i> decays at a lower order compared to the neutrino mass. We offer our 
conclusions in Sec. V. 

II. CLASSIFICATION OF EFFECTIVE AL = 2 OPERATORS 

We shall use a notation in which all fermion fields are left-handed, denoted by 



where the SU(3)c x SU{2)l x U{1)y quantum numbers are indicated in parentheses. We 
shall suppress all generation indices in this section, but will reinstate them in the next section 
when we discuss renormalizable neutrino mass models. Here L and Q stand for the lepton 
and quark doublets, respectively; = Ccr^ (C being the charge conjugation operator), 
u'^, and d'^ stand for the charge conjugates of the right-handed charged leptons, up-type 
quarks, and down-type quarks, respectively. The Standard Model Higgs doublet is denoted 
as H{1, 2, |) and H will denote its hermitian conjugate. 

The AL = 2 operators can be derived systematically using the Standard Model degrees 
of freedom as follows. There are three basic fermion bilinears that carry two units of letpon 
number: 



at least one neutrino to have a mass of about 0.03 eV |1|]. 



L(l,2 



I), e'^(l,l,l), Q(3,2,i), d%3,l,\), u%3,l,-l) 
2 b 6 6 



(1) 



(2) 
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Here i and j are SU{2)l indices, and e'^ stands for either the hermitian conjugate or the 
Dirac adjoint of e*^. Any AL = 2 effective operator will have one of these basic fermion 
bilinears accompanied by a product of other fields which is neutral under color and carries 
a net baryon number of zero. We classify the effective neutrino mass operators according 
to the number of fermion fields they contain. Three separate groups can be identified: (i) 
operators containing and no other fermion fields; (ii) operators containing four fermion 
fields; and (iii) operators containing six fermion fields. Operators containing four or more 
fermion bilinears have dimension 12 or higher and will not be considered here because, as 
mentioned above, the neutrino masses generated by such operators will be constrained by 
limits on lepton flavor violation to be typically much smaller than 0.03 eV and will not be 
that interesting for the current neutrino oscillation phenomenology. In case (i), neutrino 
masses will arise at tree level. In case (ii), one pair of fermion fields must be annihilated to 
generate neutrino masses, which will therefore arise at the one-loop level. And in case (iii), 
which requires the annihilation of two fermion pairs, neutrino masses will arise as two-loop 
radiative corrections. 

(i) With not accompanied by any more fermion fields, one obtains the well-known 
dimension five operator for neutrino mass [Q: 

O, = VUH'^H'e.keji (3) 

(ii) Operators with four fermion fields are: 

Os = Ue^W^d^H^ij (4) 



Before we list the operators in group (iii), we wish to make several remarks. 

1. Operators O2 — Oq are obtained by multiplying with any one of the combinations 
{Le'^, Q(f, Qu^i Le^, Qd^, Qu^} which all conserve lepton as well as baryon number. 
O7 and Og are obtained as the product of L*e^ with FF where F = {L, e'^, Q, u'^, d'^}. 
For operators without any derivatives or gauge boson fields (see comment 5 below), 
the number of barred (and unbarred) fields should be even. Note also that there are 
no group (ii) operators of the type e^e^ that are gauge invariant, owing to SU(2)l 
antisymmetry in the Higgs field. 

2. We have shown explicitly the SU (2) ^ group contractions using the indices ,k,l,m,n). 
The color indices are, however, suppressed. In operators involving six fermion fields 
(to follow), when there are two quark and two anti-quark fields, there are two possible 
color contractions, which will not be shown, but should be assumed. 

3. Lorentz indices are not explicitly shown in our operators. All possible Lorentz 
contractions should be allowed. Operator Oi is a unique Lorentz scalar, viz., 
{L'^^CL^)H''H^eik€ji, where C is the charge conjugation matrix. Operator O2 has 
the following Lorentz contractions: 



along with those obtained from Eq. (|^) by replacing C with cr^j,. In the first entry for 
the operator O4, the following Lorentz contractions are allowed: 



(L^^CL^)(L'=^Ce^)if'e,,efc^, iV^CL'){L^^Ce')H'e,,eM, 



(5) 



O4 = {l}^CU){Q'^Cu-)Hh,u, (V^a^MM 



-jk, 



(6) 



Although Fierz identities exist between several of these operators, they must be 
counted as independent because the SU{2)l contractions will not be the same. 
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4. We have not explicitly written down operators with AL = 2 that can be obtained by 
multiplying the lower dimensional AL = 2 operators with one of the gauge invariant, 
baryon and lepton number conserving operators that appear in the Standard Model 
Lagrangian. These latter operators include HiH^, De'^Hi, Q^dFHi, and Q'^u'^H^eij, 
as well as their hermitian conjugates. Such product operators, while not explicitly 
displayed, are understood to be present in our list. In the case of product operators 
containing HiH\ with the Higgs doublets having the trivial SU{2)l contraction as 
shown, the presence of such a higher dimensional operator will necessarily imply the 
existence of the corresponding lower dimensional operator without the HiH^ factor. 
(This can be seen by closing the Hi and W lines to form a loop, which will give an infi- 
nite contribution to the lower dimensional operator, necessitating its presence.) When 
the SU (2) i contraction on the Higgs fields is nontrivial, we do list explicitly such op- 
erators (compare, e.g., O3 with O5). Product operators involving one of the Standard 
Model operators Ue'^Hi, Q^cf^Hi, or Q^u'^H^eij (or their hermitian conjugates), and one 
of the lower dimensional AL = 2 operators of Eqs. (3)- (4), while not listed explicitly, 
can be quite interesting from the point of view of neutrino mass models. Examples 
of this type of operators are DULke''H^H'^H^eiiejm, L'UQkd''H^H"'H''euejm, and 
L^L^Q^u'^H''H"^H^eiiejjnekn- Since all possible Lorentz contractions are to be assumed, 
these operators can be generated without inducing the lower dimensional operator Oi. 
We shall give an example of this type of neutrino mass models in Sec. IIL Note 
also that the flavor structure of Ue'^Hi^ etc., in these product operators need not be 
the same as that of the corresponding Standard Model operators, which implies the 
possibility of interesting neutrino phenomenology from these operators. 

5. We have not included in our list AL = 2 operators which involve the Standard Model 
gauge boson fields. They can arise, for example, through covariant derivatives of the 
fermion or Higgs fields. It may be more difficult to generate such operators at tree 
level from an underlying renormalizable theory, making them perhaps less interesting 
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for the generation of neutrino masses. Nevertheless, we wish to hst such AL = 2 
operators of the lowest dimension, which turns out to be 7. They are {L'^ a^yL)HHB'^'' , 
{L'^a^yL)HHW>''', {L^CD^D^L)HH , and {e'=-f^D>'L)HHH . Here B^'•' and W^"" are 
the U{1)y and SU{2)l field strength tensors and we have suppressed the SU{2)l 
indices for simplicity. Although we have only shown operators with the covariant 
derivative acting on a specific field, it is understood that one should include similar 
operators with the covariant derivative acting on the other fields. For example, the 
third operator listed above includes {LFCL){D^H){D^^H), and the fourth operator 
includes {e'^Ji^L){D^H)HH , and so on. 

(iii) We now proceed to write down the operators with six fermion fields through dimen- 
sion 11. The procedure we follow is analogous to the case of the operators containing four 
fermion fields. There are 12 such operators at the dimension 9 level: 





= VL^L''e''L'e%j€ki 


Clo 




On 


= {VVQ''d''Q'd%jeku L'L^Q''d''Q'd%keji} 


Ou 




On 




Ou 


= {UUQkU'=Q''d%j, UUQiU^Q^d'eji] 




= UlJL''d%u%k 


Cl6 


= VL^e'd'e^Wsij 


On 


= VL^d^d^d^u'e.j 


Cis 


= V U d^u^u^u^Eij 


Cig 


= VQ^d'^d'^e^u^eij 


O20 


= Ud^QiU^e^u" 



And there are 40 operators with 0? = 11: 
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C»2i = {UUL^e''Q'u''H'^H^eijekmein. UU L^e''Q'u''H'^meuejmekn} 

025 = VL^Q''d''Q'u''H"'H''eim€jneki 

026 = {VVQ^d'Ue-H'H^e.iekm, VVQ^d''Lke'^H'H"'eue,^} 

027 = {VL^Q''d''Qid^H'H"'€ji€km: V Q'^d'^QJ^H' H^'eaejm} 
O2S = {VL^Q'^dVjU^H'Hieki, VL^QH'^QkU-H'H.e.i, 

VVQ''d'Qiu''H^Hiejk} 

O31 = {UL^Q,d-Q,u-H'^H'e,i, UUQ^d-Qr,u-H^H'e,ke,ie"'^] 

e»34 = e^e'^L'Q^e^d^H'^H^eikeji 
035 = e'^e'^Ue'QjU^WHh.k 
C36 = e'^e'^Q'd'^Q^d'^H'^H'eikeji 

037 = e-e-Q'd'QjU'^WHhik 

038 = e'^e'^QiU^QjU^WW 

C39 = {L^ L'' L'' LiLjH"^ H^ejm^ki: L^L^ L^L''LmLnH^H"'eijeku 

VU L^L}LiLmH^H'^ejkein, V U LpLqH^ eijekm^in^^'^} 

O40 = {VU L^Q^LiQjH"^H^ekm^ini L^L^ L'^Q''LiQiH'^H"'ejm^km 

V U l!' LiQiH^ ejm^km L^L^ L'^Q^LiQmH"^H^ejkein, 

^jk^lm 
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^ij^kl, 

C»4i = {VL^L''d%d'^H'H"'€ji€km: V U L^d%d-H' H'^eijeum} 
C»42 = {VL^L'u%u-H'H'^ejiekm, VVL''u%u-H'H'^e^,ekm} 
O43 = {VVL^d%u^H^Hiejk, L^d^LjU^H^Hieku 

VUL'^d'Uu'^H'^Hr.e^.ekme'''} 
Ou = {VL^Q''e''Qie^H'H"'€ji€km: VVQ''e''Qke^H'H'^euejm. 

VL^Q''e'Qie-H'H"'e,,ekm, VL^Q''e''Qie-H'H"'e,ke,m} 
O45 = VVe^d''e''d<'H^Hhikeji 
C46 = VL^e%''e^u^H''H^eikeji 

UUQ'Q'QkQiH'"H" 

UUQ'^Q'QkQmH'^H" 

UL^Q'^Q''QkQmH"^H"'€ii€jn, UUQ^Q^QpQqH^^H^'eijekm^in^^'^ 
UUQ^Q^'QpQqH'^H'^eikejm^ine^'^, UUQ'^Q^QpQqH'^H^^eimGjn^ki^^'^} 

048 = UUd^d'^d'^d'^H^H^eikeji 

049 = UUd''u''d'^u^H''H^eikeji 

050 = VVd''d''d^u''H^Hiejk 

052 = VL^d''u''u''u<'H''Hiejk 

053 = VVd^d%^u^HiHj 

054 = {VQ'Q'd'Q.e'^H'H^e.iekm, VQ'Q^'d'Q.e-H'H^euekm, 

VQ'Q'd'Qie'^H'H^eimejk, VQ'Q'd'Qie'^H'H^eijekm} 

055 = {VQ'QiQke^u'^H^H'eju VQ^Q^Qke-u-H^H'eiu 
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O56 = UQ^d^d^e''d<'H''H^eikeji 

058 = Vu^QjU''e''u^WH''eik 

059 = VQ^d^d^e'^u^H^Hiejk 

Oeo = Vd^QjU'^e^u'^WHi (8) 



III. RENORMALIZABLE MODELS OF NEUTRINO MASS 

The classification of tlie effective AL = 2 operators given in the previous section can 
be quite useful in building renormalizable models of neutrino mass. We shall describe in 
this section how to systematically identify from these operators interesting neutrino mass 
models. We will see that this method reproduces several well-known models. More interest- 
ingly, many new models of neutrino mass will be uncovered. While we will not present an 
exhaustive discussion of all these new models, we will outline the most interesting features 
for neutrino mass and phenomenology in several of these models. 



A. Tree level neutrino mass models 

The operator Oi of Eq. (3) can generate small neutrino masses at tree level. The 
simplest way to induce Oi is by the seesaw mechanism. As shown in Fig. 1, Oi will result 
after the heavy fields A'^i s are integrated out. Here A'^i denotes the familiar SU{2)l singlet 
right-handed neutrinos. It is also possible to induce Oi using N^, which are SU{2)l triplets 
and have zero hypercharge 

In Fig. 2, we show an alternate way of inducing Oi by the exchange of an SU (2)^ triplet 
scalar $3 which carries Y = +1. The neutral component of $3 will receive an induced 
vacuum expectation value (VEV) through its trilinear coupling with the Standard Model 
Higgs doublet. This is sometimes referred to as the Type II seesaw mechanism |jlO|, which 



11 



H H 




FIG. 1. Tree-level neutrino mass generation through Oi via the seesaw mechanism. A'^i (N^) 
is an SU{2)l singlet (triplet) fermion with zero hypercharge. 




L L 



FIG. 2. Type II seesaw mechanism where an SU{2)l triplet scalar acquires an induced VEV. 

can occur in the absence of right-handed neutrinos. For example, in S'f/(5) grand unified 
theories with a 15-plet of scalars which contain the $3 field, the requisite trilinear scalar 
coupling will arise from the Lagrangian term 15 5 5, where the 5-plet scalar fields contain 
the Standard Model Higgs doublet. The phenomenology of Oi as the source of neutrino 
mass has been studied extensively. Here we simply note that, for Oi to be interesting to the 
current neutrino oscillation data, the heavy particles Ni^^ or $3 must have masses of order 
10^2 _ iQis GeV. 
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B. One— loop neutrino mass models 



The group (ii) operators listed in Eq. (4) which involve two neutrino fields can lead to 
radiative neutrino masses at the one-loop level. Symmetries or selection rules must forbid 
the existence of the lower dimensional operator Oi for self-consistency. We shall illustrate 
five examples of one-loop neutrino mass generation. Two of the five models are well-known 
and have been thoroughly investigated in the literature, while the other three are apparently 
new. 

In this class of models, the induced neutrino masses are suppressed by a loop factor, and 
usually also by ratios of fight fermion masses to the scale of new physics A. Because of the 
additional suppression factors, the scale A in this class of models may be much smaller than 
the corresponding scale in the seesaw models, and may even be close to the electroweak 
scale. 

Notation: 

In the models we shall present in the remainder of this section, heavy scalar bosons 
will be integrated out to generate the corresponding effective AL = 2 operators. We shall 
use the following notation to denote these scalars. $ will generically denote color singlet 
scalars, fl will denote color triplet scalars, and fl color anti-triplet scalars. Note that 
is not the hermitian conjugate of Q. The SU{2)l quantum numbers of these scalar fields 
will be indicated as a subscript. Their hypercharge quantum numbers can be inferred from 
the interaction Lagrangians, and will not be displayed explicitly. For example, $i will 
transform as (1, 1) under SU{3)c x SU{2)l, while transforms as (3,3). SU{2)l indices 
will be denoted by ...), while the subscripts (a, 6, ...) will denote generation indices. 

1. Operator O2: 

Consider the following renormalizable Lagrangian: 

= fabLiLieij^^ + ii^,m2je'^ + ya4e^$2^ + h.c. (9) 
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where the Yukawa couphngs fab = —fba due to SU{2)l symmetry, ya are the (diagonal) 
Yukawa couphngs of the charged leptons, and $2 is a second scalar doublet. The hyper- 
charges of $1 and $2 are +1 and — |, respectively. The simultaneous presence of the three 
terms in Eq. (^) will result in lepton number violation, leading to the operator O21 as de- 
picted in Fig. 3. Upon closing the Le^ line (indicated by the thin solid line in Fig. 3) with 

the insertion of a Higgs field, finite neutrino masses will result. 

I 

I 



I 




FIG. 3. Diagram that induces the operator O2 as in the Zee model. 



The model just described is the Zee model for neutrino mass ITTl. The phenomenology 
of this model has been extensively studied, and we have nothing more to add here. 

2. Operator O^,'- 

The operator can be induced from the following renormalizable Lagrangian: 

= fabL\Qiei,Q, + gabVJin^i + /iOifi^^^, + h.c. (10) 

where Oi has hypercharge |, and ^2 has F = |. It is also possible to write an analogous 
Lagrangian with VLi replaced by ^3, which is a triplet of SU{2)l, rather than a singlet. 
Either of these Lagrangians will lead to via the diagram in Fig. 4. Upon closing the Q 
line with the d'^ line (indicated by the thin solid line in Fig. 4) and inserting a Higgs field, 
this diagram will induce a neutrino mass at the one-loop level. 
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H 



FIG. 4. Diagram that induces the operator O3 through the Lagrangian in Eq. dll). This 
model has a reahzation in supersymmetric models with R-parity violation. 

The induced neutrino masses can be estimated from Fig. 4 to be 

(11) 

where [mdja is the mass of the down-type quark with flavor index a, v is the VEV of H, 
and M is an average mass of the Q fields. Using f = g = 10~^, v = fi = 174 GeV and 
M = 3 TeV, we get m,^^ ~ 0.06 eV, which is in the interesting range for atmospheric neutrino 
oscillations. Assuming minimal flavor structure in the matrices / and g, we also see that 
~ {iTis/mf))'mi,.^, which is in the interesting range for solar neutrino oscillations. 
A specific realization of this model (with f^i in Fig. 4) is the supersymmetric Standard 
Model with i?-parity violation Cli is identified with d'^ and Q2 with Q in this case. 

The phenomenology of this specific realization has been well studied. We simply note that 
O3 has other realizations as well, for example, via ^3 or by using leptoquarks without 
supersymmetry. 

3. Operator O4: 

While it is possible to induce O4 by integrating out scalar fields, it turns out neutrino 
masses will arise in that case either at tree level as in Fig. 2 or at one-loop level as in the Zee 
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model depicted in Fig. 3. Consider O4 arising through the couphngs LL^-^ + Qu'^H + HH^l, 
where $3 transforms as (1, 3, +1) under SU{3)c x SU(2)l x U{1)y- In this case, the neutral 
component of $3 will acquire a tree-level VEV, and therefore a dimension 5 Majorana mass 
term for the neutrinos will be allowed at tree level, as in Fig. 2. If $3 is replaced by an 
SU{2)l singlet $1 and a second scalar doublet <l>2 is introduced so that lepton number is 
broken, then the model becomes identical to the Zee model. 

A more interesting way of generating C4 is to integrate out leptoquark gauge bosons, as 
shown in Fig. 5. The Lagrangian of the model has the form 

C-""' = Q^l,VG>t + u-l.VG^^: + GTGl^^We,, + h.c. (12) 

Here Gi is a leptoquark gauge boson which is a singlet of SU(2)l and has electric charge 
+ |. G2 is a leptoquark gauge boson that is an SU{2)l doublet with its Is = +\ member 
carrying an electric charge of +|. Upon spontaneous symmetry breaking, Gi and G2 will 
mix, leading to the generation of neutrino mass via the loop shown in Fig. 5. We have also 
allowed in Fig. 5 for the possibility that Gi may be replaced by an SU (2) l triplet gauge 
boson G3. 

I 
I 




FIG. 5. Leptoquark gauge bosons in 50(10) GUT inducing operator O4. 

In Grand Unified Theories (GUT) which unify (Q, u'^) into a common multiplet, the 
diagram of Fig. 5 will arise naturally. Such is the case in 5*0(10) GUT (but not in SU{b) 
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GUT). Consider the SU{2)l x SU{2)r x SU{A)c = G22A subgroup of 50(10). The 45 gauge 
bosons of 5*0(10) contain a (1, 1, 15) and a (2, 2, 6) under ^224- Gi of Fig. 5 is contained in 
(1, 1, 15), while G2 is in (2, 2, 6). The covariant derivative term for the 16-dimensional Higgs 
fields will contain a piece 45^ 45 16^ 16, where 45 denotes the gauge fields. The Standard 
Model Higgs doublet H may be contained partly in 16, as in currently popular S'O(IO) 
models [l^, 16 also contains a Standard Model singlet component which acquires a large 



GUT scale VEV. As a consequence, Gi and G2 will mix through this covariant derivative 
term. 

An interesting consequence of Fig. 5 is the modification of how the light neutrino masses 
scale with the quark masses in 5*0(10) models. While the conventional seesaw contributions 
to neutrino masses scale quadratically with the up-type quark masses, the new contributions 
arising from Fig. 5, given by 

--~t|^(-«).(^)_ (13) 

where g is the 50(10) coupling constant and Mq is the common mass of the leptoquark 
gauge bosons, scale linearly with the up-type quark masses. Although rriy^ here is suppressed 
by a loop factor, it may be the more dominant contribution for the lighter generations of 
neutrinos because of this linear scaling. 

4. Operator O^: 

Operator is different from 03 in that an additional HH has been added. Notice that 
the SU{2)l contraction acting on the HH is non-trivial. Here we present a renormalizable 
model which generates without inducing O3. The Lagrangian of the model is 

jC""' = fabLaQb^l + 9abLadin2 + ^2^25 H + /ifiifi^i/ + h.C. (14) 

We have not explicitly shown the SU{2)l contractions, nor the hypercharges of the colored 
scalar fields, both of which should be obvious from Eq. (0). This Lagrangian leads to 
operator through the diagram of Fig. 6. It is also possible to replace the Oi field by an 
SU{2)l triplet field fis, which is indicated in Fig. 6. 

17 



^1,3 



-A A* 
/ \ 
H H 



H 



FIG. 6. Feynman diagram that induces operator Os through Eq. (14). 



Upon closing the Qd'^ loop, we obtain a finite neutrino mass, which can be estimated to 

be 



where M is the assumed common mass of the colored scalars and {md)a stands for the mass 
of the down-type quark with flavor index a. We can obtain interesting neutrino masses from 
Eq. (15). For example, choose f = g = 10"^, A = 1, and fi = M = 5 TeV to get m^^ ~ 0.04 



5. Neutrino mass model with a product operator: 

As noted in Sec. II, operators such as Oi = UU Q^u'^H^H^-H^euejm^kni which is a 
product of the Standard Model d = A operator Q^u'^H^ekn and Oi, can lead to interesting 
neutrino mass models. Here we illustrate this possibility with Oi. The dimension five 
AL= 2 operator will not be induced in this model, since we make use of a nontrivial Lorentz 
contraction. 

Consider the following Lagrangian which breaks lepton number by two units: 




(15) 



eV. 



Oi _ 



fabLaQb^3 + gahLaUlVt2 + /ii^jfisif 



+ ^2^'2Vt'^H + ^i^Vt2^'zH + h.c. 



(16) 
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Here {Q2,^2) are SU{2)l doublet leptoquark scalars with Y = (7/6, 1/6) respectively, 
while (^3, Hg) are SU{2)l triplet leptoquarks with Y = (1/3, — 2/3) respectively. We have 
suppressed SU{2)l indices for simplicity. 

I 
I 




FIG. 7. Scalar exchange contributions inducing operator Oi (see text). 

The diagram shown in Fig. 7 will generate Oi with this Lagrangian. Upon closing the 
Q and lines, as shown by the light solid line in Fig. 7, neutrino masses will be induced 
as one-loop radiative corrections. We can estimate the induced neutrino masses to be 

where we have ignored any generation dependence in the coupling matrices / and g, and 
assumed a common mass M for the leptoquark scalars, taken to be much heavier than the 
weak scale. In Eq. (0), {mu)a denotes the mass of the up-type quark with flavor index a. 
Take, as an example, f = g = 10^^, = fi2 = fJ'S = M = 10 TeV. In this case, m^,^ ~ 0.2 
eV, which is in the interesting range for atmospheric neutrino oscillations. Clearly, other 
choices of parameters are possible with interesting neutrino phenomenology even with the 
scale of new physics being relatively low. 

Another feature worth mentioning is that the lowest dimensional AL = 2 operator that 
arises in this model is Oi as shown in Fig. 7. We can assign lepton number of —1 to both 
Q2 and ^3 that have Yukawa couplings to leptons. It is the mixing term ^22^23 that violates 
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L by two units. However, this mixing can occur in the model only after inserting three H 
fields, as in the figure. As a result, no lower dimensional AL = 2 operators are induced. 

If we replace the color anti-triplet scalars in Fig. 7, which are triplets of SU{2)l, by 
color anti-triplets that are SU{2)l singlets, the contribution to Oi will vanish, since the 
effective scalar mixing term n\n\HHH is identically zero due to SU{2) l symmetry. 



C. Two loop neutrino mass models 

Let us now turn to renormalizable models wherein neutrino masses are induced as two- 
loop radiative corrections. A generic feature of this class of models is that, for interesting 
neutrino mass phenomenology, the scale of new physics tends to be much lower than the 
corresponding scale in one-loop mass generation models. In addition to an extra loop 
suppression factor, we will see that the induced neutrino masses in this class of models scale 
quadratically with charged fermion masses, as opposed to the typical linear scaling in the 
one-loop models. 

1. Operator Ogi 



Neutrino mass models induced via Og have been well studied |jT^,|T5|. The relevant 
Lagrangian is 

C""' = fabLaL.^i + Qabelel^r + /x$i$i$r" + h.c. (18) 

where $i is a singly charged SU{2)l singlet scalar, while is a doubly charged singlet. 
Operator Og is induced through the diagram of Fig. 8. (Such a two-loop diagram for 
neutrino mass generation was first discussed in Ref. 0.) 

Upon closing the Le'^ loops on both sides, indicated by the light solid lines in Fig. 8, we 
see that the induced neutrino mass matrix has a structure given by 

Here Mi stands for the diagonal charged lepton mass matrix. An interesting feature of this 
matrix structure is that, owing to SU{2)l symmetry, fab = —fbai which means that, for 
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FIG. 8. Operator Og and the generation of two-loop neutrino masses in the model of Eq. (18). 



three generations of neutrinos, Det(mj,) = 0. This indicates that one of the neutrinos will 
remain massless at the two-loop level but it will acquire a mass at the three-loop level. 

To see the numerical magnitude of the induced neutrino masses, let us set f = g = 0.09, 
H = M = 3 TeV. The heaviest neutrino will then have a mass m^g ~ 0.03 eV. 

It should be noted that Og can also arise if $i in Fig. 8 is replaced by $3, an SU{2)l 
triplet scalar. 

2. Operator Oiq: 

We shall present two renormalizable models that induce Oiq. The effective operators are 
depicted for the two models in Figs. 9 and 10. The renormalizable Lagrangians for the two 
models can be readily written down from these figures. They are 

= fabLaU^l + QabLadin^ + KbClQ^^^ + /i^^2^2$l + k.C. (20) 

Here the two Lagrangians are meant to be taken separately, and not simultaneously. As 
shown in Figs. 9 and 10 in light solid lines, upon annihilating the L and e'^ fields as well 
as the Q and df^ fields, two-loop neutrino masses will result. The estimate of the masses is 
analogous to that in Eq. (|l^), with one difference, namely, one of the charged lepton mass 
matrices must be replaced by the down-type quark mass matrix. 
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FIG. 10. Diagram inducing Oiq through £2 ^° ■ 

It is possible to substitute an SU{2)l triplet scalar $3 for the SU{2)l singlet $1 in these 
models. This is also indicated in Figs. 9 and 10. In this case, a symmetry must be used to 
prevent the coupling ^^HH so that the neutral component of $3 does not acquire an induced 
VEV at tree level, or else the dimension 5 term of Fig. 2 will contribute more dominantly 
to neutrino masses. Clearly, both models will lead to interesting neutrino masses, e.g., with 
rriy^ ~ 0.03 eV, if the scale of new physics is around a few TeV. 

3. Operator On: 

Operator On contains the fields [LLQQd'^d'^). We shall present three models for in- 
ducing this effective operator. The three models correspond to the following different con- 
tractions: (1) {LL){QQ){d'd'), (2) {LQ){LQ){d'd^), and (3) {Ld^){Ld'){QQ). Consider the 
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{LL){QQ){d'^d'^) contraction first. The required couplings in the renormahzable model are 
shown in Fig. 11. The relevant Lagrangian for this model is 

>Cf " = fabLaU^l + QabQaQb^l + Kbdldl^^ + ^l^^U^U^ + h.C. (21) 

As shown in Fig. 11, the scalar $i may be replaced by an SU{2)l triplet scalar $3. 
In this case, the coupling ^^HH must be prevented by a symmetry so that the neutral 
component of $3 does not acquire an induced VEV at tree level, or else the dimension 5 
operator of Fig. 2 will be the more dominant effective neutrino mass operator. An example 
of such a symmetry is a Z4, under which $3 —$3, L —>■ iL, — > ie'^, with other fields 
being neutral. The cubic scalar coupling term in Eq. (l2l| ) breaks this Z4 symmetry, either 
softly or spontaneously. 

H 

I 



/ 




FIG. 11. An example of two-loop radiative generation of neutrino masses from the operator 
On. 

Upon closing the quark lines in Fig. 11, the effective scalar coupling ^^HH will be 
generated so that the neutral member of $3 will get a VEV. This VEV will be suppressed 
by two-loop factors and two powers of down-type quark masses, and will scale as v'^/M, 
where M is the assumed common mass of the heavy scalars. The neutrino masses that are 
induced will be approximately 
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m, 



) 



(22) 



(167r2)2 1^ AP 



The numerical estimate of the masses parallels that of Eq. (|19]). 



The case of using an SU{2)l singlet scalar $i, as in Eq. (|2T|), has an interesting conse- 
quence. The effective (i = 9 operator will be of the form {ve){ud){d'^d'^). In order to convert 
it into a neutrino mass term, we can contract one of the d'^ fields with the d field, but we 
will be left with {i'e){ud'^). Further exchange oi a,W boson can convert the e into a v and 
the u into a d. This is shown in Fig. 12. Note that only the longitudinal component of W 
contributes to this diagram, and we have indicated it as H~ in Fig. 12. To see this, it is con- 
venient to work in the Landau gauge, where the W boson propagator is purely transversal. 
In this gauge, it is clear that the amplitude for the transition of an off-shell vector boson 
into a scalar boson must vanish. It can, however, arise through the longitudinal component 
of the vector boson. (Such higher loop diagrams for neutrino masses have been considered 



FIG. 12. A three-loop diagram that induces neutrino masses from the operator On through 



in Ref. 0.) 



H 




H 



the Lagrangian of Eq. (^) . 



In addition to an extra loop suppression factor, the three-loop diagram of Fig. 12 
is suppressed by couplings of the longitudinal W to charged fermions. We estimate the 
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induced neutrino masses to be 



fgh ( fimjmj ' , 

where mi and rrid denote the respective charged lepton and down-type quark masses, and 
Mw is the mass of the W boson. With this estimate, we see that, for neutrino masses to be 
in the interesting range, the masses of the scalars should be of the same order as the weak 
scale. 

In Fig. 13, we display the renormalizable model where the fermion fields in On are 
contracted as {LQ){LQ){d'^d'^). The Lagrangian terms of this model are 

= fabLaQb^l + 9abdldtn[ + fiQiQiQ', + h.C. (24) 

We can also replace the field by Q3, an SU{2)l triplet. 




FIG. 13. Diagram that induces the operator On through Eq. (| 

Fig. 14 shows the fermion field contraction [Ld'^){Ld'^){QQ) that induces On. The 
Lagrangian for this case is 

= fabLadin2 + QabQaQb^l + /il^2f^2^^1 + h.C. (25) 

In both Figs. 13 and 14, the induced neutrino masses are of the same order, given 
approximately by 

- (]^) (26) 
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FIG. 14. Two-loop neutrino mass generation from On through the renormalizable Lagrangian 
inEq. 



4. Operator Ou'- 

Oyi can be obtained in a way very similar to 0\\. The neutrino mass generation mech- 
anism is shown in Fig. 15. Note that this diagram is very similar to Fig. 11, with the d"^ 
fields in Fig. 11 replaced by the u'^ fields here. We simply write down the Lagrangian for 
this model: 

i:P'' = fabLaU^l + gabQaQb^l + KbUluin, + fi^^^M + h.C. (27) 

Compared to Fig. 11, the only difference in the neutrino mass estimate is that here it will 
be proportional to m^, rather than m% 

5. Operator 02i- 

O21 is a dimension 11 operator. We shall illustrate various ways of building renormaliz- 
able models of neutrino mass with this operator. The field content of O21 is (LLLe'^Qu'^HH). 
As in the previous examples, we can contract the fermion fields in different ways. Three 
specific models are obtained by the contractions (1) {LL){e'^u'^){QL), (2) {LQ){Le'^){Lu'^), 
and (3) {LL){Lu'^){Qe^) . The Lagrangians for these three choices are as follows: 

= fabLaU^l + gabelul^i + habLaQb^l + [^i^i^iM] + h.C. 
C^^' = fabLaQb^l + 9abLael^2 + habLaUin2 + [^2^2^lHH] + h.C. 
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FIG. 15. Diagram inducing O12 through Eq. (|27|). 

Here we have not explicitly written down the renormalizable terms that result in the scalar 
mixings. These scalar mixing terms can be read off from Figs. 16, 17 and 18. 

Neutrino masses arise in all three models upon closing the Le'^ lines and the Qu'^ lines. 
The order-of-magnitude estimates for the masses in the case of Figs. 17 and 18 are 

~ i (16;^ ) [ ) 

Here we have denoted by fi all cubic scalar couplings entering the diagrams. {m^)a and {mu)a 
denote, respectively, the masses of the charged leptons and the up-type quarks with flavor 
index a. If we choose f = g = h = 0.15 and fi = M = 3 TeV, we obtain m^g ~ 0.05 eV, in 
the interesting range for atmospheric neutrino oscillations. The same estimate will apply to 
Fig. 16 if the $3 scalar is used. On the other hand, if $1 is used, there is an additional loop 
suppression, along with a suppression factor of m|/M^ as in Fig. 12. 
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FIG. 16. Generation of the d = 11 operator O21 through the Lagrangian LX^^ of Eq. (H). 
IV. NEUTRINOLESS DOUBLE BETA DECAYS AND NEUTRINO MASS 

In many of the AL = 2 effective operators, there is a tree-level contribution to neutri- 
noless double beta (/5/50i/) decay amplitudes, while neutrino masses arise only as two-loop 
radiative corrections.^ Such models will have an exciting phenomenological consequence, 
namely, that /3/30i/ decays might be observable in the current round of experiments, while 
the induced neutrino masses are quite consistent with the atmospheric and solar neutrino 
data. 

The effective operators that mediate PPOu processes will have the form [uuddee]. Oper- 
ators of this form can be easily identified from our list. The d = 9 operators from Eq. (7) 
that have this form are: 

or,l = {of!\ ofl oit 0,0}. (30) 



Gauge models where this phenomenon occurs have been noted in specific contexts |17| 
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FIG. 18. Generation of the d = 11 operator 021 through the Lagrangian 

£^21 of Eq. dH). 

Here the superscript (ii) stands for the second entry in the respective operator, and so on. 
From the set of ci = 11 operators in Eq. (8), the following operators will induce PPOu decays 
at tree level: 

w=ii _ rn(*) /o(*)'(***) n n n 

^PPOu — \^2A^ 1 ^■i2-, ^d,!^ ^■iSi 

053, oil O60}. (31) 

In order to estimate the strength of these j3j3Qv operators, we compare their amplitudes 
to that arising from light Majorana neutrino exchanges. This latter amplitude is given by 
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^/spoiy ~ G'j^rriiy (^^j, where (^^j is the average of the inverse squared Fermi momentum 
and is approximately equal to 1/(100 MeV)^. The current experimental constraint on light 
Majorana neutrino masses is m^, < 0.3 eV, derived from the nonobservation of (3(3Qv decays 
||18|| . Hence, App^i, < 10^^^ GeV^^. Consider now the amplitude arising from Oii\ Since 
this operator has dimension 9, it carries a suppression factor of A^^, where A denotes the 
scale of new physics above which the effective description becomes invalid. Comparing the 
two amplitudes, we see that, if A ~ 10^-^ GeV ~ 4 TeV, the induced PPOu decay amplitude 
will be near the current experimental limit. From the diagram of Fig. 11 (with $3 field rather 
than $1 field), we also see that rrii, ~ (l/167r^)^m^/A ~ 3 x 10~^ eV, which is extremely 
tiny and well consistent with solar and atmospheric neutrino data. (Here we considered 
only the first generation couplings and assumed the relevant / ~ (7 ~ /i ~ 1.) Similar 
estimates apply to the other d = 9 operators as well. For the d = 11 operators, there is an 
additional suppression factor of (f/A)^ due to the presence of two additional Higgs fields. 
The scale A will have to be somewhat smaller than 4 TeV in these cases for f3(30i> decays 
to be observable in the near future. The induced neutrino masses will also have this same 
additional suppression factor. Thus, it is possible to have negligible neutrino masses while 
having sizable PPOu processes. 

Of course, our estimates here have been rather crude. We have ignored important dif- 
ferences in nuclear matrix elements. However, we think that the above order-of-magnitude 
estimates suffice to demonstrate the potential importance of these operators. It will be 
interesting to investigate in more detail the rates for (3(3Qv processes in this class of models. 

V. CONCLUSIONS 

We have presented a classification of effective AL = 2 operators for the Standard Model 
that may be relevant for generating small neutrino Majorana masses. The lowest dimensional 
{d = 5) such operators are provided by the seesaw mechanism. If these d = 5 operators 
are absent from the underlying renormalizable gauge theory due to selection rules, higher 
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dimensional operators will become relevant. These higher dimensional AL = 2 operators 
will induce neutrino masses through radiative corrections. The scale of new physics, A, can 
be as low as a few TeV in this class of models. 

We have presented a list of all AL = 2 operators through d = 11 which contain no deriva- 
tives or gauge boson fields. We can readily construct from this list various renormalizable 
models for neutrino masses. We are able to identify several of the well-known radiative neu- 
trino mass models as specific realizations of some of these effective operators. Furthermore, 
we can identify a large class of new models where neutrino masses arise as one-loop, two- 
loop or even three-loop radiative corrections. We have given several examples of these new 
models, and have outlined their main features for neutrino mass phenomenology. Of special 
interest is a class of operators of dimension 9 and 11 which contribute directly to neutrinoless 
double beta decays, while generating only very tiny neutrino masses at the two-loop level. 
Models of this class have the exciting prospect that neutrinoless double beta decays will be 
observable in the ongoing round of experiments, while being fully consistent with solar and 
atmospheric neutrino oscillations. The operators presented here will have interesting conse- 
quences for cosmology as well. For example, the observed baryon asymmetry in the universe 
may have been generated by electroweak sphaleron processes which converted a primoridal 
lepton asymmetry into baryon excess. The AL= 2 operators listed here will be constrained 
if we require that the observed baryon asymmetry is produced correctly by this mechanism 
T9| . We plan to revisit our models in a future publication and analyze this and other issues 



such as the expectation for flavor changing neutral current processes in the leptonic as well 
as in the hadronic sectors. 
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